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ABSTRACT 

The Baldwin, Phillips, and Terlevich emission-line ratio diagnostic ([OIII]/H/3 versus [Nil] /Ha, 
hereafter BPT diagram) efficiently separates galaxies whose signal is dominated by star formation 
(BPT-SF) from those dominated by AGN activity (BPT-AGN). Yet this BPT diagram is limited to 
z < 0.5, the redshift at which [NII]A6584 leaves the optical spectral window. Using the Sloan Digital 
Sky Survey (SDSS), we construct a new diagnostic, or TBT diagram, that is based on rest-frame g — z 
color, [NeIII]A3869, and [OII]AA3726 + 3729 and can be used for galaxies out to z < 1.4. The TBT 
diagram identifies 98.7% of the SDSS BPT-AGN as TBT-AGN and 97% of the SDSS BPT-SF as 
TBT-SF. Furthermore, it identifies 97% of the OPTX Chandra X-ray selected AGNs as TBT-AGN. 
This is in contrast to the BPT diagram, which misidentifics 20% of X-ray selected AGNs as BPT-SF. 
We use the GOODS-N and Lockman Hole galaxy samples, with their accompanying deep Chandra 
imaging, to perform X-ray and infrared stacking analyses to further validate our TBT-AGN and 
TBT-SF selections; that is, we verify the dominance of AGN activity in the former and star formation 
activity in the latter. Finally, we address the inclusion of the majority of the BPT-comp (sources lying 
between the BPT-SF and BPT-AGN regimes) in our TBT-AGN regime. We find that the stacked 
BPT-comp source is X-ray hard ((r o ff) = LOj^o^) and has a high X-ray luminosity to total infrared 
luminosity ratio. This suggests that, on average, the X-ray signal in BPT-comp is dominated by 
obscured or low accretion rate AGN activity rather than by star formation, supporting their inclusion 
in the TBT-AGN regime. 

Subject headings: cosmology: observations — galaxies: active — galaxies: nuclei — galaxies: Seyfert 
— galaxies: distances and redshifts — X-rays: galaxies 



1. INTRODUCTION 

The most commonly used optical emission-line diag- 
nostic for separating star-forming galaxies from type II 
Active Galatic Nuclei (AGNs) relies on [OIII]A5007/H/3 
versus [NII]A6584/Ha (hereafter, BPT diagr am 
IBaldwin et alJfTMU [VeiTleux fc Osterbrockl 119871 ). The 
basic idea is that the emission lines in star-forming galax- 
ies are powered by massive stars, so there is a well- 
defined upper limit on the intensities of the collisionally 
excited lines relative to the recombination lines (such as 
Ha or H/3). In contrast, AGNs are powered by a source 
of far more energetic photons, making the collisionally 
excited lines more intense relative to the recombination 
lines. Two demarcations are commonly used for iden- 
tifying AGN-dominated gala xies versus st a r form ation 
dominated galaxies - (1) the iKewlev et all (|200lD theo- 
retical division between galaxies whose extreme ultravi- 
olet (EUV) ionizing radiation field is dominated by an 
AGN (> 50%) and those dominated by star formation 
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and (2) the iKauffmann et al.l (|2003f) empirical demarca- 
tion ba sed on the locatio n of the Sloan Digital Sky Survey 
(SPSS: [York et aLll2000f ) star-forming galaxies. Galaxies 
that lie between these two curves are often referred to 
as composite galaxies. Hereafter, we refer to these cate- 
gories as BPT-AGN, BPT-SF, and BPT-comp. 

The BPT diagram is limited in its use with optical 
spectra to galaxies with z < 0.5 (the redshift at which 
[Nil] leaves the optical spectral window). A number of 
groups have tried to extend optical emission line diagnos- 
tics to higher reds hifts by only using lines at the blue end 
of the spectrum. [Lamarcillc (2010) replaced [Nil] /Ha 
with [OII]A3726+A3729/H/3, creating the 'blue diagram', 
whi ch" can be used to classify galaxies out to z < 0.9 (see 
alsolTresse et al.lll996HRola et al.lll997tlLamareille et al.l 
120041 ^ They find that this diagnostic is very success- 
ful at identifying star-forming galaxies (with > 99% of 
the BPT-SF still classified as SF-dominated) but that it 
requires comple mentary diagnostics to robustly classify 
BPT-AGN (see lMarocco et al.|[20ll . We also note that 
[Oil] and H/3 are relatively distant in wavelength, requir- 
ing more ca reful calibr a tion a nd ex tinction corrections. 

Recently, lYan et al.l ([20111) and Uuneau et all (|2011[ ) 
introduced the CEx and MEx diagnostics, respectively, 
in which they replaced [Nil] /Ha with rest- frame U — 
B color (CEx) or with stellar mass (MEx; see also 
IWeiner et all [20 07). Both diagnostics recover the BPT- 
AGN classification very well, with 95.7% and > 99% of 
the BPT-AGN still classified as AGNs in the CEx and 
MEx diagnostics, respectively. However, because these 
diagnostics rely on H/3 at A4861 A, their use with optical 
spectra is limited to z < 1. 
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At higher redshifts the Stasihska~eFII] ((2006) DEW di- 
agn ostic, based on [4000] , [NeIII]A3869, and [Oil] (see 
also lRola et al J (19971) . can be used with optical spectra 
out to z < 1.4, at which point the lines move into the 
infrared. While these emission lines are not as strong as 
[OIII] and H/3, their ratio is the only choice for pushing 
optical spectra to these higher redshifts. [Nelll] emission 
indicates the presence of highly ionized gas and is much 
stronger than [Oil] in high-excitation AGNs. However, 
because D n [4000] requires a sufficiently high signal-to- 
noise continuum (i.e., using SDSS, we find that the uncer- 
tainty on £> n [4000] ran ges from 5% at SN/ A = 5 to 25% 
at SN/A — 1; see also lCardiel et al.lll998l) and requires 
the survey to be spectrophotometrically calibrated, its 
usage with distant galaxies is limited. 

In this article we examine whether rest-frame g — z 
color, which requires near-infrared imaging at the higher 
redshifts, is a compelling replacement for £) n [4000]. 
While the ratio of [NeIII]/[OII] alone effectively separates 
metal-rich star-forming galaxies from AGNs, metal-poor 
star- forming galaxies have high values of [NeIII]/[OII] (as 
a result of less line blanketing enabling a harder stellar ra- 
diation field). Fortunate ly, metal-poor galaxie s also tend 
to be bluer (see Fig. 7 in lTremonti et al.l[200 4), so we can 
use their color to distinguish them from AG Ns (which 
tend to be bulge dominated and redder, see lYan et all 
120111 ). We refer to this new diagnostic as the TBT dia- 
gram. 

We first establish the reliability of the TBT di- 
agnostic in reproducing the BPT classifications at 
low redshifts using SDSS, the largest spectroscopic 
sample to date of emission-line galaxies. We then 
test the TBT diagnostic at higher redshifts using 
the highly spectroscopicall y complete O PTX X-ray 
selected sample of AG Ns (ITrouille et ail 120081 [20091: 
iTrouille fc Bargerl 120101 ) . The misidentification of X- 
ray selected AGNs as SF-dominated galaxies is a po- 
tential issue for all optical emi s sion-li ne diagnostic di- 
agrams. In ITrouille fe B"argen (pOlCl ) we found that 
~ 20% of the Lx > 10 42 erg s _1 Chandra X-ray selected 
AGNs in our OPTX sample that have [OIII], H/3, [Nil], 
and Ha fluxes with signal-to-noise greater than five are 
misidentified by th e BPT diagram as BPT-SF (see also 
iWinter et alll2010l for evidence of this in the Swift BAT 
sample and iBongiorno et al.ll2010l for evidence of this in 
the XMM-COSMOS sample). Existing high-redshift op- 
tical emission-line diagnostics also misidentify a signif- 
icant fraction of X-ray selected AGNs, e.g., ~ 8% and 
- 22% of X-ray selected AGNs in the MEx and CEx 
diagrams lie in the MEx-SF and CEx-SF regimes, re- 
spectively. 

Stacking techniques have been widely used in X-ray 
astronomy to study the average properties of source 
populations selected to have certain well-defined prop- 
erties and which are too X-ray fain t to be detected 
indiv i dually (e.g.. iBrandt et all 120011: lAlexander et all 
20011: iHornschemeier et all l2002t iNandra et al.l 120021: 



Georgakakis et all 20031: Laird et al.l 120051 : iLehmer et all 
20051 l2008t iTreister et all 120091) . As a proof-of-concept 



of the TBT diagnostic, we perform an X-ray stacking 
analysis of the TBT-SF and TBT-AGN in the highly 
spectroscopically complete Great Observatories Origins 
Deep Survey North (GOODS-N: IGiavalisco et"all[200l 



with accompany i ng Ch andra Deep Field North (CDFN; 
lAlexander et al.l (2003) imaging to determine whether 
they are, on average, X-ray soft or X-ray hard. 

An X-ray hard source is indicative of obscured AGN 
activity or the presence of high-mass X-ray binaries 
(HMXBs) associated with ongoing star formation. An 
X-ray soft source is indicative of unobscured AGN activ- 
ity or the presence of low-mass X-ray binaries (LMXBs) 
associated with old stellar populations. In order to 
distinguish between these scenarios, we perform an in- 
frared (IR) stacking analysis using the Spitzer Space 
Telescope 24/im data. Numerous studies have found 
a relation betwen the X-ray and IR luminosities for 
star-forming galax i es and for AGNs dPtak et all 120031: 
Persic et all 120041: [Alexander et all 120051: iTeng et al l 
20051 : IGeorgakakis et all 120071) . We compare the X-ray 
and IR luminosities for our stacked TBT-SF and stacked 
TBT-AGN with the results from these studies to verify 
the dominance of star formation activity in the former 
and AGN activity in the latter. 

We then perform X-ray and IR stacking analyses of 
the BPT categories to confirm the presence of AGN 
activity in BPT-comp, as suggested by the TBT diag- 
nostic. In most optical emission-line diagnostics, the 
BPT-comp lie in a transition regi on or within the SF- 
dominated regime. For example, iKewlev et al.l (2006) 
find that BPT-comp are indistinguishable from HII re- 
gions and are significantly distinct from Seyferts in an 
[OIII] /[Oil] versus [OI] /Ha plot. In both the 'blue di- 
agram' and the CEx diagnostic, the majority of BPT- 
comp lie within their SF-dominated regimes (~ 83% and 
~ 75%, respectively). In contrast, the bulk of BPT-comp 
lie within our TBT-AGN regime (see also the MEx di- 
agnostic). Because BPT-comp are a significant percent- 
age of the overall low-redshift emission-line galaxy pop- 
ulation (e.g., in SDSS3, BPT-SF, BPT-comp, and BPT- 
AGN comprise 69%, 20%, and 11%, respectively, of the 
overall population), inclusion or exclusion of BPT-comp 
in AGN samples can have an important impact on re- 
sults. In order to have the necessary statistics required to 
do a robust stacking analysis, we combine our GOODS- 
N galaxy sample with galaxy samples from two Lockman 
Hole (LH) fields, all of which have deep Chandra imag- 
ing. 

The structure of the paper is as follows. In Section [2] 
we briefly describe the SDSS sample, our OPTX X-ray 
selected sample of AGNs, and our optical spectroscopic 
samples of the GOODS-N/LH fields. In Section 1 we 
use the SDSS sample to calibrate our TBT diagnostic 
to match the BPT diagram. In Section 0] we determine 
how well our TBT-AGN selection matches with an X- 
ray selection of AGNs. In Section [5] we do X-ray and 
IR stacking analyses to verify the reliability of our TBT- 
AGN and TBT-SF classes. In Section © we do X-ray 
and IR stacking analyses on our BPT-comp to determine 
whether the implications of the TBT diagnostic with re- 
spect to the dominance of AGN activity in BPT-comp are 
confirmed. In Section [7] we compare our TBT diagnostic 
to other diagnostics for separating star-forming galaxies 
from AGNs. In Section [5] we summarize our results. 



6 Here we have restricted the SDSS DR8 sample to emission-line 
galaxies whose [Nil], Ha, [OIII], and H/3 fluxes have a signal-to- 
noise ratio greater than five. 
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All magnitudes are in the AB magnitude system. We 
assume Ojvf = 0.3,51a = 0.7, and Ho = 70 km s _1 
Mpc^ 1 . 

2. SAMPLE 

2.1. SDSS: Low-Redshift Galaxy Sample 

SDSS has obtained deep, multi-color images covering 
more than a quarter of the sky with follow-up spec- 
troscopy of over a million objects. Here we use the 
SPSS spectroscop ic data from Data Release 8 (DR8; 
lAihara et al.ll20Tlh . We use the emission-line fluxes mea- 
sured by the MPA-JHU group as described in Section 
4.3 of the data release paper. These fluxes are esti- 
mated from simultaneous Gaussian fits to the continuum 
subtracted spectra to account for stellar absorption and 
line blending. DR8 provides spectra for 868,492 different 
galaxies. The SDSS spectral range is 3800 - 9200 A. 
Since the [Oil] line lies at 3726 A and the [Nil] line 
lies at 6584 A, we limit our SDSS samples to galaxies 
with 0.02 < z < 0.35. Our SDSS BPT sample consists 
of the 243,865 SDSS galaxies that have Ha, [Nil], H/3, 
and [OIII] fluxes with signal-to-noise ratio (SNR) greater 
than five (out of the 818,333 spectra for different galaxies 
in the DR8 sample). Of these SDSS BPT galaxies, 23,048 
also have both [Nelll] and [Oil] fluxes with SNR > 5. We 
refer to these as our SDSS TBT sample. 

2.2. OPTX: X-ray Selected Sample of AGNs 

The OPTX sample consists of 1789 X-ray sources in 
two intermediate depth wide-field surveys in the Lock- 
man Hole region of low galactic c olumn density and 
one deep pencil-beam survey (CDFN; iBrandt et al.ll200il : 
lAlexander et al.l |20 03 ). The Lockman Hole fields are the 
Chandra Large A rea Synoptic X-ray Survey (CLASXS; 
lYang et al. 20QH) and the Chandra Lockman Area North 
Surve y (CLANS; iTrouille et ai1[2008l [20091 : IWilkes et al.l 
2009). We have spectroscopically observed 84% of the 
OPTX sources using th e DEep Imaging Multi-Object 
Spectrograph fDEIMOS; fFaber et al.lf2003D on the Keck 
II 10 m telescope and the HYDRA multi-object spec- 
trograph on the WIYN 3.5 m telescope (for details of 
the o bservations and reduction process, see lTrouille et al.l 
l2008|) 



In ITrouille et a l. (2008. 120091) we used the X-ray fluxes 
and spectroscopic redshifts to calculate rest-frame 2 — 
8 keV luminosities, Lx- At z < 3 (which is all we 
consider in this paper), we calculated the luminosities 
from the observed-frame 2 — 8 keV fluxes, assuming 
an intrinsic X-ray spectral index of F = 1.8. That 
is, Lx = / x &xd\ x fc— correction, where for z < 3, 
fc— correction = (1 + z) T ~ 2 and / = fu-s keV- Using the 
individual X-ray spectral indice s (e.g., (r) ~ 1.75 with a 
dispersion of ss 0.33 derived bv lTozzi et all 120061 for the 
X-ray bright CDFS sources), rather than the universal X- 
ray spectral index of F = 1.8 adopted here, to calculate 
the A;— corrections would result in only a small difference 
(an average factor of 0.9) in the rest-frame luminosities. 
We have not corrected the X-ray luminosities for absorp- 
tion since these corrections a re small in the 2 — 8 keV 
band (e.g.. lBarger et al.1l2002D . and we are only using the 
X-ray luminosities to identify sources as X-ray AGNs. 

In the following, we limit our study to the 561 OPTX 
X-ray sources with spectroscopic redshifts, whose 2 — 



8 keV flux has a significance greater than 3 er, and whose 
Lx > 10 42 erg s _1 (247, 163, and 151 sources from 
the CLANS, CLASXS, and CDFN fields, respectively). 
Lx > 10 42 erg s" 1 is a co mmonly used conservative 
threshold for AGN activity (iHornschemeier et al.l 120011: 
Barger et aHl2002t iSzokolv et al.l 120041: iSilverman et al l 



2005 1 ICoiletalE oOQ) that is based o n energetic grounds 
dZezas et al.lll998tlMoran et al.lfl999h . Using the calibra- 
tion bv iRanalli et al.1 ((2003), one would need a star for- 
mation rate (SFR) of 200 M Q yr _1 to produce enough 
X-ray l uminosity from non-AGN s to cross this thresh- 
old. In ITrouille fc Barged (|2010( ) we found a < 5% con- 
tamination rate of our OPTX sample by sources with 
SFR > 200 Mp, yr" 1 . We determined this using the 
iMagnelli et ail pOOll space densities at a range of red- 
shifts for ultraluminous infrared galaxies (ULIRGs) with 
estimated SFR > 172 M Q . 

Our study focuses on optical emission-line ratio diag- 
nostics based on flux ratios. We comp ute the relative 
line fl uxes from the spectra using the iTremonti et al.l 
(2004) software. In brief, we subtract the stellar con- 
tinuum and absorption lines by fitting a linear combina- 
tion of single stellar popu lation models of different ages 
(jBruzual fc Charlotl l200l . We remove any remaining 
residuals from the continuum using a sliding 250 A me- 
dian. The relative line fluxes and errors are estimated 
from simultaneous Gaussian fits to the continuum sub- 
tracted spectra. 

Our DEIMOS spectra are of high quality from ~ 
4800 A to - 9300 A, such that the [Oil] and [Nelll] lines 
lie within our spectral window from z ~ 0.3 — 1.4. Of our 
561 OPTX X-ray sources, 197 are non-broad-line AGNs 
(FWHM < 2000 kms" 1 ; hereafter non-BLAGNs) and lie 
within this redshift range. We only include non-BLAGNs 
in this analysis, since the narrow lines in BLAGNs are 
overwhelmed by the emission from the broad-line region. 
We identify 103 OPTX X-ray selected AGNs that have 
[Nelll] and [Oil] fluxes with SNR > 5. 

2.3. GOODS-N/LH: Higher- Redshift Galaxy Sample 

Each of the three OPTX fields is the focus of a com- 
prehensive spectroscopic follow-up of all galaxies in the 
field. The CDFN encompasses the intensively studied 
GOODS-N field, which we describe in more detail below. 
The CLANS and CLASXS fields reside in the Lockman 
Hole (LH) region of low Galactic column density. 

A random spectroscopic sample of 3082 K S ,A& < 21.5 
galaxies in these two LH fields was obtained using 
DEIMOS on Keck II (L. Cowie, priv. comm.). We here- 
after refer to this combined sample as the LH galaxy 
sample. 

The GOODS-N field has among the deepest images 
ever obtained in a number of bandpasses, including the 
Chandra 2 Ms CDFN image. It also has been the tar- 
get o f extensive spectros c opic observations over the years 
200(1 iWirth et all 120041: iCowie et al.l 
2008; ICooper et al.l 120111 ). With its 



,e.g., iCohen et al.l 
2004; H - 
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high optical spectroscopic completeness and deep X-ray 
coverage, the GOODS-N catalog provides an ideal sam- 
ple for studying the average X-ray properties of optically 
s elected samples of A GNs . 

Bar ger et al.l (|2008l ) presented a highly complete spec- 
troscopic survey of the GOODS-N field. Over the years a 
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number of groups have made observations of this region, 
first primarily us ing the Low-Res olution Imaging Spec- 
trograph (LRIS; lOke et al.l 119951 ) on the Keck I 10 m 
telescope and later using DEIMOS. Barger et al. added 
to these samples by observing all missing or unidenti- 
fied galaxies to date with DEIMOS. In order to pro- 
vide a uniform spectral database, they also re-observed 
sources where the original spectra were of poor quality 
or where previous redshifts were obtained with instru- 
ments other than DEIMOS. The redshift identifications 
are now greater than 90% complete to magnitudes of 
F435W-AB = 24.5 and K StAB = 21.5 and to 24 /im fluxes 
of 250 /tJy. The final Barger et al. catalog provides spec- 
troscopic redshifts for 2710 galaxies in this field. 

We compute the relative line fluxes for the relevant 
emission lines in the combined GOODS-N/LH galaxy 
sample using the same method and software as de- 
scribed for the OPTX X-ray sample above. To cre- 
ate our GOODS-N/LH BPT diagram, we use the 727 
sources that have [OIII], H/3, [Nil], and Ha fluxes with 
SNR > 5. For our GOODS-N TBT diagram, we use 
the 670 sources that have [Nelll] and [Oil] fluxes with 
SNR > 5. As discussed in the previous subsection, be- 
cause of the DEIMOS spectral window, the GOODS-N 
sources plotted in our TBT diagnostic have a redshift 
range from z = 0.3 — 1.4. 

In Figure Q] we show the redshift distributions (out to 
z = 1.5) for the SDSS, GOODS-N, LH, and OPTX sam- 
ples. The large SDSS sample has a median redshift of 
(z) ~ 0.1, whereas the GOODS-N, LH, and OPTX sur- 
veys allow us to test our TBT diagnostic out to z = 1.4. 
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Fig. 1. — Spectroscopic redshift distribution out to z = 1.5 for 
(a) the SDSS BPT sample, (b) the GOODS-N galaxy survey, (c) 
the LH galaxy survey, and (d) the OPTX survey. The different 
redshift ranges seen in (a)-(c) reflect the magnitude limits of the 
different surveys. 



tometr y into g and z magnitudes at z = using kcorrect 
v4_l_4 (|Blanton fc Roweisl [2007). For z = 1.4 sources, 
the g— band (~ 0.5/im) redshifts into the J— band (~ 
1.2/tm) and the z— band (~ 0.9/tm) redshifts into the 
K— band (~ 2.2/mi). As a result, fields with ex- 
tensive photometric coverage use observed-frame band- 
passes that are very close to the rest-frame g— and z— 
bandpasses of interest, and so the fc— corrections are 
small. For our TBT diagnostic, we subtract the rest- 
frame z— band from the rest-frame g— band to obtain the 
rest-frame g — z color; hereafter, (<? — z). 

To determine the error associated with the (g — z) 
color for each sample, we randomly alter the magnitudes 
by an amount consiste nt with the photometric uncer- 
tainties and re-run the IBlanton fc Roweisl (120071) soft- 
ware. We then determine the 1 a standard deviation 
on A = (g - ^original - (.9 - z) ra ndom- For our 



SDSS sample, a A = 0.21. For our GOODS-N/LH sam- 
ple, <ta = 0.07. The high-quality photometry and exten- 
sive coverage of our GOODS-N/LH fields ensures these 
lo w fc— correction err ors. 

iChilingarian et all ([2010) compare the 

Blan ton fc Roweisl (120071 ) kcorrect code with their 
own to explore the systematic error associated with 
fc— correction software. For fc— corrected optical colors 
(e.g., [g — r], [r — z]), the differences between the 
resulting colors from the two codes are relatively small 
(~ 0.05 or less). Therefore, for our study, systematic 
errors are not significant. However, we note here for 
completeness that for fc— corrected colors in which one 
band is in the rest-frame IR (e.g., [r - H], [r - K]), 
the differences can be significant (~ 0.15; see their 
Fig. 7), and the systematic error need be taken into 
account. 

3. A NEW Z < 1.4 EMISSION-LINE RATIO DIAGNOSTIC 
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2.4. Rest-frame g, z Magnitudes 

We have u,q,r.i,z m agnitudes for the SDSS galaxies 
([Abazajian et al.l l2009 ! ). g, r, i, z, J, 

H,K S magnitudes for the OPTX AGNs an d LH 
galaxies (|Trouille et al.l l2008t iKeenan et al.ll2010D . and 
B,V,RJ,Z \ J,H,K„ magni t udes for the GOODS- 
N galaxies (iGiavalisco et all I2004t iBarger et all |2008t 
IKeenan et al.l I2010IV We transform the observed pho- 



FlG. 2.— BPT diagram - log([NII]/Ha) vs. log( [OIII] /H/3) - 
for the SDSS galaxies that have [OIII], H/3, [Nil], and Ha with 
SNR > 5. The dashed curve shows the Kauff mann et al.l l|2003) 
empirical division betwe en star-forming gala xies and AGNs. The 
dotted curve shows the Kcwlcy et al. (2001) theoretical division. 
As discussed in the text, we represent the SDSS data using a com- 
bination of a two-dimensional histogram and plotted points. BPT- 
SF are shown in blue, BPT-comp in gray, and BPT-AGN in red. 

Figure [2] shows the BPT diagram (based on flux ra- 
tios of the specified lines) applied to our SDSS BPT 
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Fig. 3.— (a) TBT diagram - (g-z) color vs. log([NcIII]/[OII]) 
- for the SDSS galaxies that have [OIII], H/3, [Nil], Ha, [Nelll], 
and [Oil] fluxes with SNR > 5. The top, middle, and bottom pan- 
els show the positions of the BPT-AGN, BPT-comp, and BPT- 
SF, respectively. The dashed line provides the empirical separa- 
tion maximizing the fraction of BPT-AGN to total population in 
the upper-right of the diagram (see Eq. [TJ. (b) Same as (a) but 
for SDSS galaxies that have [OIII], H/3, [Nil], Ha, and [Oil] with 
SNR > 5, but [Nelll] has SNR < 5. For the 15% with [Nelll] < 0, 
we set [Nelll] equal to the 1 a error. 



sampl e. The dotted curve indicates the iKewlev et al.l 
(120011) theoretical divi s ion an d the dashed curve indicates 
the IKauffmann et al.l (|2003l ) empirical division between 
AGNs and star-forming galaxies, as discussed in the In- 
troduction. In this and subsequent figures we represent 
the SDSS data using a combination of a two-dimensional 
histogram and plotted points. We histogram the data 
where more than four data points fall in an individual 
pixel (the size of the pixel is determined by dividing the 
plot into 150x150 bins) and plot it as individual points 
otherwise. The histograms have been square root scaled 
for better visibility. Each step in contour level represents 
20% fewer sources. The darkest contours correspond to 
391 sources, 42 sources, and 151 sources for the BPT-SF, 
BPT-AGN, and BPT-comp, respectively. 

We find that 69% of our SP SS BPT sample lie be- 
low the IKauffmann et al.l (120031) division in th e BPT-SF 
regime, 11% lie above the IKewlev et al.l (|2001[ ) division, 
in the BPT-AGN regime, and 20% lie in between the two 
divisions, in the BPT-comp regime. If instead we ap- 
ply the BPT diagnostic to our SDSS TBT sample (i.e., 
SDSS galaxies that have [OIII], H/3, [Nil], Ha, [Nelll], 
and [Oil] fluxes with SNR > 5), we find that 60% are 
BPT-SF, 32% are BPT-AGN, and 8% are BPT-comp. 

The BPT diagram can only be used to identify AGNs 
out to z = 0.5, the redshift at which [Nil] leaves the 
optical spectral window. The ratio of [Nelll] to [Oil] is 
a good candidate for pushing optical narrow emission- 
line ratio diagnostics to higher redshifts because both 
lines are relatively strong, lie in the blue end of the spec- 
trum (measurable in optical spectra out to z < 1.4), and 
are close in wavelength (thus avoiding reddening effect 
issues). Furthermore, [Nelll] has a significantly higher 
ionization potential (63.45 eV) than [Oil] (35.12 eV). As 
a result, the [NeIII]/[OII] ratio is higher in AGNs than in 
star-forming galaxies. We also note that the ionization 
potential for [Nelll] is significantly higher than for [OIII] , 
whose ionization potential is 54.94 eV. In Section [7] we 
discuss how AGN selection using a [NeIII]-based diag- 
nostic may be more discerning than one based on [OIII] 
as a result of this higher ionization potential (requiring 
a harder ionizing flux). 

In Figure [3ja) we plot our TBT diagnostic - (g — z) 
color versus log([NeIII]/[OII]) - for our SDSS TBT sam- 
ple, color-coded according to the BPT classifications in 
Figure O The top, middle, and bottom panels show 
the locations of the BPT-AGN, BPT-comp, and BPT- 
SF, respectively. The darkest contours correspond to 71 
sources, 16 sources, and 68 sources for the BPT-AGN, 
BPT-comp, and BPT-SF, respectively. The BPT-SF 
show a trend with color in that the bluer BPT-SF exhibit 
higher values of [NeIII]/[OII]. This provides a separation 
in color-space between the BPT-AGN and the BPT-SF 
with high values of [NeIII]/[OII]. This trend with color 
is likely due to the fac t that bluer galaxie s are more 
metal poor (see Fig. 7 in lTremonti et al. 2004) and hence 
have harder ste llar radiation fields (higher [NeIII]/[OII]; 
iShi et al.ll2007l) . On the other hand, few BPT-AGN re- 
side in ver y blue galaxies. A s discussed in detail in Sec- 
tion 3.3 of lYan et all (|201lD . nearly all BPT-AGN are 
found in red galaxies or in galaxies with intermediate 
colors between red and blue. 
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The dashed line designates 



o.o 



(g-z) = -1.2 x log([NcIII]/[OII]) - 0.4 , (1) 



our empirical separation maximizing the fraction of 
BPT-AGN to total population in the upper-right of the 
diagram. Hereafter, we refer to the sources in the upper- 
right (lower-left) of our TBT diagnostic as TBT-AGN 
(TBT-SF). 

We find that 98.7% of the BPT-AGN lie in the TBT- 
AGN regime and 97.2% of the BPT-SF lie in the TBT-SF 
regime. Likewise, we find that only 3.5% of the sources in 
the TBT-AGN regime are BPT-SF and 1% of the sources 
in the TBT-SF regime are BPT-AGN. The BPT-comp 
lie on either side of the division, with 68.6% in the TBT- 
AGN regime. Overall the BPT-comp constitute 5% and 
16% of the TBT-SF and TBT-AGN, respectively. 

There are 200,712 galaxies in our SDSS BPT sam- 
ple that have [Oil] with SNR > 5 but [Nelll] with 
SNR < 5. In Figure [3(b) we show the TBT diagram 
for these sources. For the 15% with [Nelll] < 0, we set 
[Nelll] equal to the 1 a error. The darkest contours cor- 
respond to 94 sources, 170 sources, and 546 sources for 
the BPT-AGN, BPT-comp, and BPT-SF, respectively. 
We find that ~ 85% of the BPT-SF lie within the TBT- 
SF regime, ~ 86% of the BPT-AGN lie within the TBT- 
AGN regime, and ~ 62% of the BPT-comp lie within the 
TBT-AGN regime. The trends follow those for our SDSS 
TBT sample. The 8% of SDSS BPT galaxies for which 
neither [Oil] nor [Nelll] has SNR > 5 and the < 1% of 
SDSS BPT galaxies for which [Nelll] has SNR > 5 but 
[Oil] has SNR < 5 are not considered here. 

We have trained our TBT diagnostic on the SDSS 
galaxies, which have (z) ~ 0.1. Because our TBT di- 
agnostic can be used with optical spectra out to z = 1.4, 
we examine the impact of metallicity evolution with red- 
shift on our empirical separation between TBT-AGN 
and TBT-SF. Galaxy metallicities decrease by a fac- 
tor of ~ 0.3 dex between th e local val ue and the value 
at z ~ 2 (lErb et al.l 120061: see also ICowie fe Bargerl 
12001 IKewlev fc Ellison! I2008t IZahid et al.l 1201 ID . For 
an M star = 10 10 M Q galaxy, this corresponds to a 
shift from 12 + log(0/H) = 8.6 to 8.3 (note that 
more massive galaxies undergo less metallicity evolu- 
tion). Using the lShi et al.l (|2007T ) relation between metal- 
licity and log([NeIII]/[OII]), we find that, in this case, 
log([NeIII]/[OII]) shifts by only - 20%, from -0.89 to 
—0.72. Applying this 20% increase in the ratio of [Nelll] 
to [Oil] to all SDSS BPT-SF, we find that an additional 
5% move to the TBT-AGN regime. This corresponds to 
an increase of only 5.5% in the number of TBT-AGN that 
are BPT-SF. Similarly, there is only a 1.8% increase in 
the number of TBT-AGN that are BPT-comp. Because 
the impact is relatively small, in this article we do not 
consider any metallicity evolution with redshift in our 
empirical separation between TBT-SF and TBT-AGN. 

We also considered the impact of color evolution on our 
TBT diagnostic. At higher redshifts, galaxies are bluer 
as a result of higher speci fic star formation rates. P urely 
passive evolution models (iBruzual fc Cha riot 2003) with 
an instantaneous burst and a lChabrierl (|2003f) initial mass 
function predict a A (g — z) ~ 0.24 between z = and 
1.4, for a formation redshift of 5. Applying this color evo- 
lution to all SDSS BPT-SF, the impact is in our favor. 



The BPT-SF move down the y-axis in our TBT diagnos- 
tic to lie even further below our empirical separation, in 
the TBT-SF regime. Applying this color evolution to all 
SDSS BPT-AGN, we find that an additional 3.6% move 
to the TBT-SF regime. This corresponds to an increase 
of 2.7% in the number of TBT-SF that are BPT-AGN. 
Similarly, there is only a 2.9% increase in the number 
of TBT-SF that are BPT-comp. Because the impact is 
small, in this article we do not consider any color evo- 
lution with redshift in our empirical separation between 
TBT-SF and TBT-AGN. 



4. 



COMPARING THE TBT DIAGNOSTIC WITH AN X-RAY 
SELECTION OF AGNS 
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Fig. 4.— (a) TBT diagram - (g — z) color vs. log([NeIII]/[OII]) 
— for the X-ray selected non-BLAGNs in our OPTX sample with 
0.3 < Z < 1.4 for which the [Nelll] and [Oil] fluxes have a 
SNR > 5. Small black triangles — non-BLAGNs with 10 42 < 
Lx < 10 44 erg s _1 . Large black triangles — non-BLAGNs with 
Lx > 10 44 erg s" 1 . The dashed line provides the empirical sepa- 
ration maximizing the fraction of BPT-AGN to total population in 
the upper-right of the diagram (see Figure l3l and Eq.[TJ. (b) Same 
as (a) but for the OPTX X-ray selected non-BLAGNS for which 
[Oil] has a SNR > 5 but [Nelll] has a SNR < 5. For the 6% with 
[Nelll] < 0, we set [Nelll] equal to the 1 a error. 



In lTrouille fc Barged poToj) we found that only a little 
over half (52%) of the X-ray selected non-BLAGNs in our 
OPTX sample lie in the BPT-AGN regime of the BPT 
diagram. These sources form a sequence similar to that 
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of the BPT-AGN, emerging from the HII region sequence 
and extending to the upper-right hand side of the BPT 
diagram. 

Of the X-ray selected non-BLAGNs in our OPTX 
sample, 20% are misi dentified as BPT-SF, i.e. , as pure 
star- forming galax ies (jTrouille fc Bargcr 120101 see also 
iWinter et alfelOl for evidence of this in the Swift BAT 
sample). I ncreased extinction does n ot account for these. 
Instead, in iTrouille fc B argcr (2010) we note that these 
misidentified sources have lower i[oiii]/^Jf ratios than 
those that lie in the BPT-AGN regime. We postulate 
that the low forbidden emission line strengths in the 
misidentified sources are a result of the complexity of the 
structure of the narrow-line region, which causes ionizing 
photons from the central engine to not be absorbed. 

The misidentiheation of X-ray selected AGNs as star- 
forming galaxies is a potential issue for all optical 
emission-line diagnostic diagrams, including our TBT 
diagram. In Figure IHa) we plot the TBT diagram - 
(5 - z) color versus log([NeIII]/[OII]) - for the 0.3 < 
z < 1.4 X-ray selected non-BLAGNs in our OPTX sam- 
ple that have [Nelll] and [Oil] with SNR > 5. We find 
that 97% (100/103) of our X-ray selected AGNs lie in 
the TBT-AGN regime. Thus, the TBT diagnostic does 
a much better job of correctly identifying X-ray selected 
AGNs than the BPT diagnostic, misidentifying only 3% 
as TBT-SF compared to the BPT diagnostic's misiden- 
tiheation of 20% as BPT-SF. 

As mentioned in Section [2~2l there are 94 0.3 < z < 1.4 
OPTX non-BLAGNs that do not fulfill the criteria of 
having both [Nelll] and [Oil] with SNR > 5. Sixteen 
of these have neither [Oil] nor [Nelll] with a SNR > 5. 
In Figure BJb) we show the remaining 78 OPTX X-ray 
selected non-BLAGNs for which [Oil] has a SNR > 5 but 
[Nelll] has a SNR < 5. No sources have only [Nelll] with 
a SNR > 5. We hnd that 92% (72/78) lie in the TBT- 
AGN regime. This supports our results for the OPTX 
sources in which both [Nelll] and [Oil] have a SNR > 5. 

5. VERIFICATION OF THE TBT DIAGNOSTIC: STACKING 
ANALYSES 

The GOODS-N sample of spectroscopically observed 
galaxies with deep Chandra and Spitzer imaging out to 
z = 1.4 provides an ideal pilot for conhrming the dom- 
inance of AGN activity in the TBT-AGN and star for- 
mation activity in the TBT-SF. 

Figure [5] shows our TBT diagnostic for the 642 0.3 < 
Z < 1.4 GOODS-N sources for which the [Nelll] and [Oil] 
fluxes have a SNR > 5. We have excluded sources with 
broadened emission lines (FWHM > 2000 km s _1 ) from 
this analysis. The dashed line indicates our empirical 
division between the TBT-AGN and TBT-SF (see Figure 
Hand Eq.rjJ. We And that 189 GOODS-N sources lie in 
the TBT-AGN regime, corresponding to a space density 
of 4.6 x 10~ 5 Mpc~ 3 optically selected GOODS-N TBT- 
AGN. 

There are 11 GOODS-N sources in our TBT diagram 
that are cross-listed as non-BLAGNs in the CDFN cata- 
log with Lx > 10 42 erg s _1 (black triangles). As dis- 
cussed in Section [4] our TBT diagnostic has a much 
lower rate of misidentifying X-ray selected AGNs as star- 
forming galaxies than the BPT diagnostic (~ 3% versus 
20%). We use our OPTX X-ray selected sample of AGNs 
for that analysis. Since the GOODS-N field lies within 
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Fig. 5.— TBT diagram - (g - z) color vs. log([NeIH]/[OH]) - 
for the GOODS-N 0.3 < 2 < 1.4 sources for which the [Nelll] and 
[Oil] fluxes have a SNR > 5. Large black triangles — non-BLAGNs 
cross-listed in the Chandra catalog with Lx > 10 42 erg s . Or- 
ange diamonds — non-BLAGNs cross-listed in the Chandra catalog, 
but with Lx < 10 42 erg s — 1 . These sources are significantly de- 
tected in both the 0.5 — 2 keV band and 2 — 8 kcV bands. Magenta 
diamonds — non-BLAGNs cross-listed in the Chandra catalog, but 
with L x < 10 42 erg s _1 . These sources are significantly detected 
in the 0.5 — 2 keV band but not the 2 — 8 keV band. Black cir- 
cles — remaining sources. The dashed line provides the empirical 
separation maximizing the fraction of BPT-AGN to total popula- 
tion in the upper-right of the diagram (see Figure [3] and Eq. [TJ. 



the CDFN, one of our OPTX fields, it is no surprise that 
only one of the L x > 10 42 erg s -1 non-BLAGNs (black 
triangles) in Figure [S] is misidentified as a TBT-SF. 

In the following section we perform X-ray and IR stack- 
ing analyses in order to investigate whether the remain- 
ing L x < 10 42 erg s" 1 TBT-AGN harbor obscured AGN 
activity, verifying the validity of our TBT-AGN selec- 
tion. We also perform stacking analyses of our TBT-SF 
in order to verify the dominance of star formation activ- 
ity in these sources. The effective X-ray spectral slope 
(r c ff ; derived from the ratio of the 2 — 8 keV band to the 
0.5 — 2 keV band) is indicative of different source types. 
Unobscured AGNs exhibit T ff > 1.7. In obscured AGNs, 
photoelectric absorption of X-ray soft photons by neutral 
gas along the line of sight causes a flattening of the slo pe, 
such that r cff < 1.7 (although see lTrouille et al.ll2009l for 
a discussion of discrepancies between X-ray and optical 
spectral type). At the same time, X-ray binaries exhibit 
a range in X-ray spectral slopes. LMXBs, which are as- 
sociated with old stellar populations, tend to be softer in 
X-rays, with P > 1.7. HMXBs, which are associated with 
ongoing star formation, tend to be hard er in X-rays, with 
T = 0.5 - 1 dColbert et al.ll200l. See iFabbianol (|2006f ) 
and lRemillard fc McClintockJ (|2006h for detailed reviews 
of X-ray binary populations. In order to distinguish be- 
tween these scenarios, we perform an IR stacking analysis 
in Section I5~2l 

We restrict our analysis to lower redshift sources (z < 
0.7) in order to not be misled by the automatic soften- 
ing of the effective X-ray spectral slope as one observes 
to higher redshifts (i.e., as the redshift increases, the 
0.5 — 2 keV and 2 — 8 keV bands are sampling higher ener- 
gies whose photons can more easily penetrate obscuring 
material) . 
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5.1. X-ray Stacking 

Quantitatively, if stacking a sample of 'n' objects yields 
'JV,' counts in a signal aperture of area l A s ' and l iV&' 
smoothed background counts in that same area l A s \ then 
the mean number of source counts per object in the signal 
region is 

(N) src =-[N S -N b \. (2) 
n 

We obtain an estimate of the mean spectral slope of 
the detected signal by performing the stacking in two 
energy bands, soft (0.5 — 2 keV) and hard (2 — 8 keV), 
and deriving an effective power-law photon index, (r o g). 

Here we use the S TACKFAST X-ray stacking program 
(jHickoxet al.|[2007t) . In STACKFAST, 'A s ' is denned as 
the area enclosed within rgo from the source position, 
where rgo is an approximation of the 90% point-spread 
function (PSF) energy encircled radius at 1.5 keV, and 
varies atQ: 

r 90 = 1" + 1O"(0/1O') 2 , (3) 

with 9 equal to the off-axis angle. In order to maximize 
the number of source counts, rath er than limit i ng th e 
stacking to only the central 6' as in IHickox et all (|2007h . 
we use the central 10' around the pointing center for each 
observation. 

A few bright sources would dominate our estimate 
of the mean spectral shape, so we exclude from our 
analysis sources that lie close to or are associated with 
a known X-ray detected source, hereafter call e d co- 
incidental contaminan ts (see also IHickox et al.l 120071 : 
iGeorgantopoulos et alJl2008t iFiore et all 120081 ). To this 
end, we first applied a mask to all known X-ray detected 
source positions. We used 3 x rgo (see Eq. 3) as the mask 
aperture radius. 

We cross reference the l A s ' for our sources of inter- 
est with the X-ray photon locations from each Chandra 
pointing (i.e., each OBSID event file). We then create 
smoothed 0.5 — 2 keV and 2 — 8 keV background maps 
using the CIAO WAVDETECT task and determine the 
background counts within the same t A s \ Subtracting 
this background and dividing by the number of sources 
being stacked, we obtain the average X-ray signal in 
counts per source. We create exposure maps using cus- 
tom routines (A. Vikhlinin, private communication) and 
determine the total exposure time for each source being 
stacked. Dividing the total counts in the stacked source 
by the total exposure time, we derive the average count 
rate (counts per second) for the stacked source. We use 
the ratio of the 2 — 8 keV count rate to the 0.5 — 2 keV 
count rate to derive (r c ff). 

Count uncertainties are calculated using the approxi- 
mation \/ X + 0.7 5+ 1, where X is the number of counts 
m a given band ()Gehrelslfl986l) . We set a significance 
threshold of 3 a. Uncertainties in the hardness ratio and 
(r c ff ) are derived by propagating these count rate errors. 

Table [T] shows the results from our stacking analysis 
for the GOODS-N TBT-SF and TBT-AGN. Column 2 
lists the total number of sources in each of these cate- 
gories. Column 3 provides the total number of sources 
used in the stacking analysis, after excluding individually 
X-ray detected sources and coincidental contaminants. 

1 http:/ /cxc. harvard.edu/proposer/POG 



Columns 4 and 5 state fluxes and detection significance 
for the given X-ray band. We consider that stacked 
sources with < 3 a detection are not significantly de- 
tected. Column 6 provides the (T e g) value for the stacked 
source. In column 7 we list the volume- weighted redshift 
for the stacked source and in column 8 we list the de- 
rived 2 — 8 keV luminosity based on the /2-8 keV an d the 
volume- weighted redshift. 

We provide specifics for each category in the following 
subsections. 

5.1.1. TBT-SF 

There are 155 z < 0.7 GOODS-N TBT-SF. As dis- 
cussed above, we restrict our analysis to these lower red- 
shifts in order to not be misled by the automatic soften- 
ing of the effective X-ray spectral slope as one observes 
to higher redshift. Only one of these z < 0.7 GOODS-N 
TBT-SF is directly X-ray detected, and then only in the 
0.5-2 keV band. Using the 2 - 8 keV flux limit for the 
CDFN image (see Table [2J , we determine an upper limit 
to its hardness ratio and find that it is X-ray soft, with 

r > i.7. 

After excluding this known X-ray source, as well as co- 
incidental contaminants, we use the STACKFAST pro- 
gram to determine the average X-ray signal in the re- 
maining 148 TBT-SF. The stacked source is significantly 
detected in the 0.5 — 2 keV band (4.7 a) but not in the 
2 - 8 keV band (2.2 a). We find a (r off ) = 1.5+g;J, al- 
though we caution that this is based on a low-significance 
2 — 8 keV signal. Given the large uncertainties, the 
stacked source could be X-ray soft (r > 1.7). 

5.1.2. TBT-AGN 

There are 72 z < 0.7 GOODS-N TBT-AGN. Six are 
directly X-ray detected in both the 0.5 — 2 keV and 2 — 
8 keV bands. Three of these are obvious X-ray selected 
AGNs with Lx > 10 42 erg s _1 . The three remaining X- 
ray detected sources have Lx < 10 42 erg s _1 . All three 
are X-ray hard, with T < 1.4. 

After excluding these known X-ray sources, as well 
as coincidental contaminants, we use the STACKFAST 
program to determine the average X-ray signal in the 
remaining 54 TBT-AGN. The stacked source is signifi- 
cantly detected in both bands (5.4 a in the 0.5 — 2 keV 
band and 4.0 a in the 2 - 8 keV band). We find 
a (r o ff) = l.Olgg, consistent with being X-ray hard 
(r<1.4). 

5.1.3. Monte-Carlo Simulation 

We carried out a series of Monte Carlo (MC) stacking 
simulations to assess false-detection probabilities empir- 
ically. For each category we performed 1000 trials and 
used the same number of stacked sources and the same 
procedure as in the original stacking, albeit with random 
RA and Dec positions. 

In Section [5TTT1 we found that the TBT-SF are signif- 
icantly detected in the 0.5 — 2 keV band but not in the 
2-8 keV band. In Section [5XU we found that the TBT- 
AGN are significantly detected in both the 0.5 — 2 keV 
and 2 — 8 keV bands. 

Our MC simulations yield a 0.1% probability of gen- 
erating the observed 0.5 — 2 keV flux for our stacked 
TBT-SF source and a 52% probability of generating the 
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TABLE 1 

TBT X-ray Stacking Analysis Results 



Category 
(1) 


Total # 
(2) 


# used in stack 
(3) 


r a 
J 0.5-2 keV 

(4) 


f a 
/2-8 keV 

(5) 


<reff> 
(6) 


z b 
(7) 


logZ/2-8 kov [erg s L ] 
(8) 


TBT-SF 


155 


148 


0.62 ±0.13 (4.7 a) 


2.50 ± 1.1 (2.2 a) 


1 c+0.7d 


0.58 


40.53 d 


TBT-AGN C 


72 


54 


1.13 ± 0.2 (5.4 a) 


5.87 ± 1.4 (4 a) 


1 +0 ' 3 


0.58 


40.86 



2 In units of 10 



-17 



,-2 .-1 



°Valumc-wcighted redshift for the stacked source. 
c Only including sources with z < 0.7. 
We caution that this is based on a low significar 



keV signal. 



SFR (Mg/^year) 




44.0 



44.5 



45.0 
LOG L(TIR) 



45.5 



46.0 



Fig. 6. — 1/2-8 kev/^TiR versus log Ltir for our 0.3 < z < 
0.7 GOODS-N galaxies that have [Nelll] and [Oil] fluxes with 
SNR > 5. Red filled squares - directly X-ray detected TBT-AGN. 
Large (small) symbols indicate sources with Lx > 10 42 erg s — 1 
(Lx < 10 42 erg s" 1 ). Blue filled star — directly X-ray detected 
TBT-SF. The downward pointing arrow indicates that the source 
is undetected in the 2 — 8 keV band and is assigned the 2 — 8 keV 
flux limit (see Table [2} • The upward and leftward pointing arrows 
indicate that the source is undetected in the Spitzer 24/im image 
and is assigned the 24/^m flux limit (see Table |2J. Red open square 
- stacked TBT-AGN, excluding directly X-ray detected sources. 
Blue open star — stacked TBT-SF, excluding directly X-ray de- 
tected sources. Solid line - expected ratio for the HMXB popula- 
tion in a ga l axy w ith the given SFR (see top axis), according to 
I Persic ct al. (2004); dashed line - expected ratio for the overall star 
formation occurrin g in a galaxy with the given SFR, according to 
IMineo et afl (t20TTD . 

observed 2 — 8 keV flux. Our MC analysis confirms that 
the TBT-SF clearly have an excess of 0.5 — 2 keV counts, 
well above those we obtain randomly. The 2 — 8 keV sig- 
nal is within the noise. Thus, the observed X-ray softness 
in our TBT-SF appears to be a reliable result. 

Our MC simulations for the TBT-AGN yield a 0% 
probability that we would see the observed 0.5 — 2 keV 
flux or the observed 2 — 8 keV flux. Our MC analysis 
confirms we are recovering a real signal in both bands. 
None of our MC simulations for the TBT-AGN result in 
both the 0.5 — 2 keV and 2 — 8 keV bands detected at 
a 3 a level or greater. In other words, there is a ~ 0% 
false-alarm rate for our stacked TBT-AGN signal. Thus, 
the observed X-ray hardness in our TBT-AGN appears 
to be a reliable result. 

5.2. IR Stacking 

In Section [5. II we derived the average X-ray properties 
for the GOODS-N TBT-SF and TBT-AGN using an X- 
ray stacking analysis. The stacked TBT-AGN source was 
significantly detected in both the 0.5—2 keV and 2—8 keV 
bands and was quite X-ray hard, with (T c g) = LOlgJ. 



TABLE 2 
Flux Limits 



Chandra 2 — 8 keV Image 


Flux Limit ( erg cm '* s L ) 


CDFN 


1.5 x 10-"' 


CLANS 


3.5 x 10" 15 


CLASXS 


3.5 x 10~ 15 


Spitzer 24/im Image 


Flux Limit (/^Jy) 


GOODS-N 


75 


LH a 


150 



"Includes both the CLANS and CLASXS fields. 

This X-ray hard signal could be due either to AGN activ- 
ity with obscuration or HMXBs associated with ongoing 
star formation. 

Numerous studies have found that in star-forming 
galaxies without an AGN, the t otal X-ray luminos- 
ity correlat e s wi t h the SFR dNandra et al.l 120021: 



Baue r et al 
20031: 



Hornschemeie 



120021: iRanalli et all 
Colbert et al.l 12004 



Rovilos et al.l 



et 



2T 



2009). 



2003: Grimm 



_ iPersic et al.. 

.20051: IPersic fc Rephaeli 
Furthermore. iDavid et al 



et al.1 



2004; 



2007; 
(1992) 

found a linear relation between the far-infrared (FIR) 
luminosity and the 0.5 — 4.5 keV luminosities for a 
sampl e of starburst galaxies ob ser ved by Einste i n (see 
also iFabbiano fc Shaplevi I2002D . IRanalli et all (|2003f ) 
extended this study to the 2—10 keV band using 
the A SCA and BeppoSAX satellites. Persi c et al.l 
(2004) determined the relation between the X-ray 
and FIR luminosities for both the full contribution 
from star formation activity and for HMXBs alone. A 
number of groups have shown that AGNs and AGN- 
starburst composi t e galaxies lie ab o ve th e se relations 
dPtak et all [200l I Alexander et al.l [20051: iTeng et al.l 
I2005t iGeorgakakis et al.l 120071) . The additional X-ray 
luminosity is attributed to AGN activity in these 
galaxies. 

Following these previous studies, here we use the IR 
properties of our TBT-AGN to determine whether the 
X-ray hard signal is due to AGN activity with obscu- 
raiton or HMXBs. We also check whether the X-ray 
soft signal in our T BT-SF is cons i stent with pure star- 
forming galaxies. In lTrouille et al.1 ((20091) we determined 
the Spitzer 24/im fluxes and luminosities for our OPTX 
X-ray selected AGNs. We follow the same procedure 
here to determine /24/im and £24^111 f° r our GOODS-N 
galaxies. 

To determine /24/im for our stacked TBT-AGN and 
stacked TBT-SF, we use th e publicly available IAS Stack- 
ing Library IDL software (|Bethermin et al.ll201Q[ ). The 
IAS Stacking software uses a DAOPHOT-type photom- 
etry IDL procedure, APER, with a preset PSF for the 
Spitzer 24/zm band (13" for the object aperture and sky 
radii of 20" — 32"). A median stacking is preferable to 
mean stacking because the median analysis is more sta- 
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ble and robust to small numbers of bright sources. Using 
the volume-weighted redshift for our stacke d source, we 
transf orm /24^ m into L24^m- We then use the lRieke et al.l 
(2009) eq. A6 to transform L24jum m t° the total infrared 
luminosity, L^ir- 

The IAS Stacking software provides the 1 a standard 
deviation on the stacked flux. However, given that the 
Spitzer image resolution is low, we need to consider the 
likelihood of misidentifications and overlap. To assess 
how well this 1 a standard deviation reflects the contam- 
ination from spurious signals, we carry out a series of 
Monte Carlo stacking. We perform 1000 trials and use 
the same number of stacked sources and follow the same 
procedure as in the original stacking. The only difference 
is that we use random RA and Dec positions. For both 
our stacked TBT-AGN and our stacked TBT-SF, > 99% 
of our simulations result in stacked fluxes less than our 
1 a error. While < 10% of the random RA, Dec positions 
in each simulation do overlap with real 24/xm sources (as 
expected given the low image resolution), because we use 
a median stacking, these spuriously high fluxes are ex- 
cluded from the stacked signal. 

In Figure [6] we plot the ratio of the X-ray to total IR 
luminosities, Lx/Ltir, versus I/tir for our TBT cate- 
gories. The top axis shows t h e asso ciated SFR for a given 
Zvtir, following iRieke et al.l (I2009T) . The solid lin e shows 
the expected ratio derived by iPersic eFall (|200l for the 
HMXB population in a galaxy with the gi ven SFR. The 
dashed line shows the iMineo et al.l (j 2 1 lh expected ra- 
tio for all star formation activity in a galaxy with the 
given SFR. The stacked TBT-SF (blue open star) was 
significantly detected in the 0.5 — 2 keV band (4.7 cr) 
but not in the 2 — 8 keV band (2.2 a). Nonetheless, 
we use the derived Lx from the stacking analysis (see 
Table [T]) and determine the 1 a error on Lx / £tir by 
propogating the errors on both the stacked X-ray and 
IR signal. The stacked TBT-SF Lx/Ltir is consistent 
with that expected for pure star-forming galaxies. This 
corroborates what we find for the one individually X- 
ray detected TBT-SF (blue filled star). Since this source 
is undetected in both the Spitzer 24/im image and the 
CDFN 2 — 8 keV image, we assign it the flux limits for 
these images (see Table [2]) and use the arrows to des- 
ignate it as corresponding to upper limits. Given these 
uncertainties, its Lx/Ltjr is not inconsistent with pure 
star-forming galaxies. 

The three X-ray detected TBT-AGN with L x > 
10 42 erg s _1 (large red filled squares) and the three X-ray 
detected TBT-AGN with L x < 10 42 erg s" 1 (small red 
filled squares) lie clearly above the expected range for 
Lx/Ltir for pure star-forming galaxies. The stacked 
TBT-AGN (red open square) also lies well above this 
range. We determine the 1 a error on Lx/^tir f° r the 
stacked source by propagating the errors on both the 
stacked X-ray and IR signal. We find that the stacked 
TBT-AGN lies > 3 a above the expected range for 
Lx/Ltir for pure star-forming galaxies, supporting our 
hypothesis that, on average, TBT-AGN harbor AGN ac- 
tivity. 

6. BPT DIAGNOSTIC: STACKING ANALYSES 

In Section [3] we found that the majority of the SDSS 
BPT-comp lie within the TBT-AGN regime. Here our 
primary goal is to do X-ray and IR stacking analyses to 
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Fig. 7.— (a) BPT diagram for the GOODS-N/LH z < 0.5 sources 
that have [OIII], H/3, [Nil], and Ha fluxes with SNR > 5. Large 
black triangles — non-BLAGNs cross-listed in the Chandra cata- 
log with Lx > 10 42 erg s —1 . Orange diamonds — non-BLAGNs 
cross-listed in the Chandra catalog, but with Lx < 10 42 erg s _1 . 
These sources are significantly detected in both the 0.5 — 2 keV 
band and 2 — 8 keV bands. Magenta diamonds — non-BLAGNs 
cross-listed in the Chandra catalog, but with Lx < 10 42 erg s _1 . 
These sources are significantly detected in the 0.5 — 2 keV band 
but not the 2 — 8 keV b and. Black circles — rema ining sources. The 
dashed curve shows the Kauffmann et al. (2003) empirical division 
bet ween star-forming ga laxies and AGNs. The dotted curve shows 
thejKcwlcy ct al. (200T]) theoretical division. 

investigate whether the BPT-comp signal, on average, is 
dominated by AGN or star formation activity. 

Since the BPT diagram is restricted to sources with 
z < 0.5, there are too few GOODS-N BPT-comp for a 
robust stacking analysis (specifically, there are only 22 
GOODS-N BPT-comp). To increase our numbers, here 
we also include the LH galaxy sample. Figure[7]shows the 
BPT diagram for the 727 z < 0.5 GOODS-N/LH galaxies 
that have [OIII], H/3, [Nil], and Ha fluxes with SNR > 5. 
We have excluded sources with broadened emission lines 
(FWHM > 2000 km s" 1 ) from this analysis. 

As discussed in Section^ the BPT diagnostic does not 
match 100% with an X-ray selection of AGNs. In Fig- 
ure [7] we see that 50% (9/18) of the X-ray selected non- 
BLAGNs with Lx > 10 42 erg s _1 (black triangles) are 
BPT-AGN, 28% (5/18) are BPT-comp, and 22% (4/18) 
are BPT-SF. 

In the following sections we do stacking analyses to 
verify that the remaining Lx < 10 42 erg s _1 BPT- 
AGN and BPT-SF are, on average, AGN-dominated and 
SF-dominated, respectively. Having tested our stacking 
analyses on these more secure BPT categories, we apply 
them to our BPT-comp. 

6.1. X-ray Stacking 

Table [3] provides the results from our X-ray stacking 
analysis for the GOODS-N/LH BPT-SF, BPT-AGN, and 
BP T-com p. We followed the same procedure as in Sec- 
tion 15.11 We provide specifics for each category in the 
following subsections. 

6.1.1. BPT-SF 

There are 605 GOODS-N/LH BPT-SF. Twenty-one 
are directly X-ray detected. Four of these are obvious 
X-ray selected AGNs with L x > 10 42 erg s -1 . Of the 17 
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TABLE 3 

BPT X-ray Stacking Analysis Results 



Category 
(1) 


Total # 
(2) 


# used in stack 
(3) 


/0.5 - 2 keV a 
(4) 


/2-8 koV a 

(5) 


(Teff) 

(6) 


z b 
(7) 


logZ/2-8 kcv [erg s 1 ] 
(8) 


BPT-SF 


605 


448 


0.66 ±0.12 (5.5 a) 


2.82 ± 1.23 (2.3 a) 


1 K+U.7C 

■■■^-Q.a 
0.72l^ c 


0.41 


40.20 c 


BPT-AGN d 


31 


17 


5.39 ±3.37 (1.6 a) 


39.4 ± 22.9 (1.8 a) 


0.41 


41.35 c 


BPT-Comp 


91 


76 


1.71 ±0.33 (5.2 a) 


8.92 ± 2.23 (4.0 a) 


1 0+ ' 4 

■■■• U -Q.4 


0.41 


40.70 



1 111 



.its of 10' 



17 , 



"Volume- weighted redshift for t 

Wc caution that this is based on a low significa 
d Three out of four GOODS-N BPT-AGN and 11 

Stacking the remaining 17 GOODS-N/LH X- 
e Wc caution that this is based on a low signif 



ay 



s 2 - 8 kcV signal, 
t of 27 LH BPT-AGN a 
ndetected BPT-AGN pn 
3 0.5 - 2 kcV and 2-8 



detected : 
>t statistics 



for the sta 



X-ray detected sources with L x < 10 42 erg s _1 , ten are 
only detected in the 0.5—2 keV band. Using the 2 — 8 keV 
flux limit for their Chandra image (see Tabled, we de- 
termine an upper limit to their hardness ratios and find 
that all ten are X-ray soft, with T > 1.7. Of the seven 
remaining sources that are detected in both bands, three 
are X-ray soft, with T > 1.7. Therefore, the majority 
(13/17) of the X-ray detected, Lx < 10 42 erg s _1 sources 
are X-ray soft, with T > 1.7. 

After excluding these known X-ray sources, as well as 
coincidental contaminants, we use the STACKFAST pro- 
gram to determine the average X-ray signal in the re- 
maining 448 BPT-SF. The stacked source is significantly 
detected in the 0.5 — 2 keV band (5.5 a) but not in the 
2 - 8 keV band (2.3 a). We find a (T) = 1.5±j^, al- 
though we caution that this is based on a low-significance 
2 — 8 keV signal. Given the large uncertainties, the 
stacked source could be X-ray soft (r > 1.7). 

Our MC simulations (see Section [5 .1.31 for details of the 
procedure) yield a 2% probability of randomly generating 
the observed 0.5 - 2 keV flux for our stacked BPT-SF 
and a 96% probability of generating the observed 2 — 
8 keV flux. Our MC analysis confirms that these BPT- 
SF clearly have an excess of 0.5 — 2 keV counts, well 
above those we obtain randomly. The 2 — 8 keV signal is 
within the noise. 

6.1.2. BPT-AGN 

There are 31 GOODS-N/LH BPT-AGN. Fourteen are 
directly X-ray detected. Nine of these are obvious X-ray 
selected AGNs with L x > 10 42 erg s" 1 . Of the five X- 
ray detected sources with Lx < 10 erg s _1 , one is only 
detected in the 0.5 — 2 keV band. Using the 2 — 8 kcV flux 
limit for its Chandra image (see Table [2]) , we determine 
an upper limit to its hardness ratio and find that it is X- 
ray soft, with T > 1.7. The four remaining sources that 
are detected in both bands are X-ray hard, with T < 1.4. 

After excluding these known X-ray sources, as well 
as coincidental contaminants, we use the STACKFAST 
program to determine the average X-ray signal in the 
remaining 17 BPT-AGN. Given the small number of 
sources in this stacking analysis, it is not surprising that 
the stacked source is neither significantly detected in the 
0.5-2 keV band (1.6 a) nor in the 2-8 keV band (1.8 a). 
We find a (r) = 0.72^02, although we caution that this 
is based on a low-significance 0.5 — 2 keV and 2 — 8 keV 
signal. Given the large uncertainties, this stacked source 
could be X-ray hard (r < 1.4). 

Our MC simulations (see Section 15.1.31 for details of 
the procedure) yield a ~ 98% probability of randomly 



generating the observed 0.5 — 2 keV and 2 — 8 keV fluxes. 
Our MC analysis confirms that the 0.5 — 2 keV and 2 — 
8 keV signals are within the noise. 

6.1.3. BPT-comp 

There are 91 GOODS-N/LH BPT-comp. Eleven are 
directly X-ray detected. Five of these are obvious X-ray 
selected AGNs with Lx > 10 42 erg s _1 . Of the six X-ray 
detected sources with Lx < 10 42 erg s _1 , three are only 
detected in the 0.5 — 2 keV band. Using the 2 — 8 keV flux 
limit for their Chandra image (see Table[2]), we determine 
an upper limit to their hardness ratios and find that all 
three are X-ray soft, with T > 1.7. The remaining three 
sources are detected in both bands and are X-ray hard, 
with r < 1.4. 

After excluding known X-ray sources, as well as coinci- 
dental contaminants, we use the STACKFAST program 
to determine the average X-ray signal in the remaining 
76 BPT-comp. The stacked source is significantly de- 
tected in both bands (5.2 a in the 0.5 — 2 keV band and 
4.0 a in the 2 - 8 keV band). We find a (T eS ) = 1.0t° A 4 , 
i.e., the source is X-ray hard. 

Our MC simulations for the BPT-comp yield a 5% 
probability that we would see the observed 0.5 — 2 keV 
flux and a 2% probability we would see the observed 
2 — 8 keV flux, confirming that we are recovering a real 
signal in both bands. Fewer than 1% of our MC simula- 
tions for the BPT-comp result in both the 0.5 — 2 keV and 
2—8 keV bands detected at a 3 a level or greater. In other 
words, there is a < 1% false-alarm rate for our stacked 
BPT-comp signal. Thus, the observed X-ray hardness in 
our BPT-comp appears to be a reliable result. 

6.2. IR Stacking 

As discussed in Section 15.21 the X-ray hard signal in 
our stacked BPT-comp could be due either to HMXBs 
or to AGN activity with obscuration. We follow the 
same procedure as in Section 15.21 to distinguish between 
these scenarios using the IR. We determine the 1 a error 
on Lx I £tir for the stacked BPT-comp (gray open dia- 
mond) by propagating the errors on both the stacked X- 
ray and IR signal. Figure |S] shows that the stacked BPT- 
comp lies > 3 a above the expected range for Lx / £tir 
for pure star- forming galaxies (dashed line). This corrob- 
orates what we find for the individually X-ray detected 
BPT-comp (gray filled diamonds), where the majority 
(8/11) have Lx/Ltir values well above the expected 
range for pure star-forming galaxies. These results sup- 
port our TBT diagnostic inclusion of the bulk of BPT- 
comp in the TBT-AGN regime. 
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We also note that every one of the BPT-AGN, BPT- 
SF, and BPT-comp with L x > 10 42 erg s" 1 (large filled 
symbols) have Lx / £tir values well above the expected 
range for pure star-forming galaxies. This further con- 
firms that the signal in these X-ray selected AGNs is 
dominated by AGN activity. An optical classification as 
BPT-SF is inaccurate. 

The stacked BPT-AGN (red open square in Figure [8]) 
was not signitificantly detected in cither the 0.5—2 keV or 
2 — 8 keV bands (1.6 a and 1.8 a, respectively). Nonethe- 
less, we use the derived Lx from the stacking analysis 
(see Table and determine the 1 a error on Lx / £tir 
by propagating the errors on both the stacked X-ray and 
IR signal. The stacked BPT-AGN L x /L Tm lies well 
above the expected range for pure star-forming galaxies. 
This corroborates what we find for the individually X-ray 
detected BPT-AGN with L x < 10 42 erg s" 1 (small red 
filled squares), where the majority (3/5) have Lx/Ltir 
well above the expected range for pure star-forming 
galaxies. The two with Ltir > 10 44 5 lie within the 
expected range for pure star-forming galaxies. These 
two sources are an additional example of the BPT di- 
agnostic potentially misidentifying sources - in this case, 
SF-dominated sources as BPT-AGN. 

The stacked BPT-SF (blue open star in Figure [5]) was 
significantly detected in the 0.5 — 2 keV band (5.5 a) 
but not in the 2 — 8 keV band (2.3 a). Nonetheless, 
we use the derived Lx from the stacking analysis (see 
Table |3|) and determine the 1 a error on L x / Ltir by 
propagating the errors on both the stacked X-ray and 
IR signal. The stacked BPT-SF Lx/L^m is consistent 
with that expected for pure star-forming galaxies. This 
corroborates what we find for the individually X-ray de- 
tected BPT-SF with L x < 10 42 erg s" 1 (small blue filled 
stars), where the majority (12/17) have Lx/^tir values 
consistent with the expected range for pure star-forming 
galaxies. 



SFR (Ma/year) 



7. DISCUSSION 



COMPARISON WITH ALTERNATIVE 
DIAGNOSTICS 



As mentioned in the Introduction, lLamareillel ()2010f ) 
investigate the use of a 'blue diagnostic' - [OIIIJ/H/3 ver- 
sus [OII]/H/3 - to classify galaxies as s tar-forming galax- 
ies, AGNs, and composites. Although lLamareillel (|2010l ) 
separate out LINERa3, here we include both AGNs and 
LINERs in the BPT-AGN category, as in the rest of our 
article. Since [OIII] , [Oil] , and H/3 lie at the blue end of 
the spectrum, this diagnostic can be applied to galaxies 
with optical spectra out to z ~ 0.9. This method pro- 
vides a successful means for creating a reliable sample 
of star-forming galaxies out to these redshifts. Their SF- 
dominated regime (SFG in their Table 1) encloses > 99% 
of the BPT-SF. Only 1.5% of the sources in the SFG 
regime are BPT-AGN. The 'blue diagnostic', however, 
is not as effective in creating a reliable AGN selection. 
Their combined Sy2, SF/Sy2, and LINER regimes iden- 
tify - 94% of all the BPT-AGN. However, ~ 23% of the 

7 The right wing of the BPT diagram further subdivides into 
an upper and lower branch, with LINERs being t he sources in the 
lower branch exhibitin g lower [OIII] luminosities ( Hcckmanl ll9801 : 
Kauffmann ct al. 2003). The nature of these sources (weak AGNs 
versus 'retired' galaxies dominated by old stellar populations with 
relatively hard radiation fie l ds) is still much deba ted (Kcwle y et al.l 
2006; Stasihska ct al. 2008; Cid Fcrnandcs ct al. 2010, 2011). 
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Fig. 8. — L2—8 kev/LTiR versus log Ltir for our z < 0.5 
GOODS-N/LH galaxies that have [OIII], H/3, [Nil], and Ha fluxes 
with SNR > 5. Red filled squares - directly X-ray detected 
BPT-AGN. Blue filled stars - directly X-ray detected BPT-SF. 
Gray filled diamonds — directly X-ray detected BPT-comp. Large 
(small) symbols designate sources with Lx > 10 42 erg s _1 (Lx < 
10 42 erg s — 1 ). Downward pointing arrows indicate sources unde- 
tected in the 2 — 8 keV band (but detected in the 0.5 — 2 keV band) 
that are assigned the 2 — 8 ke V flux limit (see Table [2]l . Upward 
and leftward pointing arrows indicate sources undetected in the 
Spitzer 24fim image that are assigned the 2Afim flux limit (see Ta- 
ble[2j). Blue open star — stacked BPT-SF, excluding directly X-ray 
detected sources. Red open square - stacked BPT-AGN, exclud- 
ing directly X-ray detected sources. Gray open diamond - stacked 
BpT-comp, excluding directly X-ray detected sources. Solid line 
— expected ratio for the HMXB pop ulation in a galaxy with the 
given SFR (see top axis), according to Persic ct al. (2004). Dashed 
line - expected ratio for the overall star formation occurrin g in a 
galaxy with the given SFR, according to Minco et al. (20fj|). 

sources in this combined regime are BPT-SlfJ 

A classification sc heme based on [NeIII]/[OII] 
( Stasihska et al.l 120061 : this work) complements this 
'blu e diagnost i c' by correctly classifyi ng galaxies in 
the lLamareillel ()2010D SF/Sy2 category. iMarocco et~aT| 
( 20111) recently examined the loc ation of the | Lamareillel 
(2010) SF/Sy2 galaxies in the iStasihska et al.l (|2006l ) 
DEW diagnostic - [NeIII]/[OII] versus D„[4000]. They 
find that by applying an additional empirical cut in 
the DEW diagnostic space to these SF/Sy2 galaxies, 
they are able to correctly identify 99% of the BPT-SF 
and 97% of the Seyfert 2s. Overall, this combined 
diagnostic correctly identifies 85% of Seyfert 2s and 99% 
of BPT-SF. 

Our TBT diagnostic, based on [NeIII]/[OII] versus 
rest- frame g — z color (rather than D n [4000] ) , results in 
minimal overlap between the BPT-SF and BPT-AGN, 
with only ~ 1.3% of the BPT-AGN lying within the 
TBT-SF regime and only 2.8% of BPT-SF lying within 
the TBT-AGN regime. Likewise, of the sources in the 
TBT-AGN regime, only 3.5% are BPT-SF and of the 
source s in the TBT-SF reg ime, only 1% are BPT-AGN. 

The lJuneau et all ([2TJT1I ) MEx diagnostic - [OIII] /H/3 
versus stellar mass - also complements the 'blue diagnos- 
tic' in its ability to reliably identify AGNs. The MEx- 
AGN regime encloses 99% of BPT-AGN. Only 6% of the 

8 If we remove the LINER regime from this analysis, we note 
that the Sy2 and SF/Sy2 regimes enclose 59% and 26% of the 
Seyferts, respectively. While < 3% of sources in the Sy2 regime 
are BPT-SF, 74% of sources in the SF/Sy2 regime are BPT-SF. 
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sou rces in the MEx- AGN regime are BPT-SF. Similarly, 
the lYan et ail ([20111 ) CEx diagnostic - [OIII]/H/3 versus 
rest-frame U — B color - also complements the 'blue di- 
agnostic', with the CEx-AGN regime enclosing 95.7% of 
BPT-AGN. Only 1.9% of the sources in the CEx-AGN 
regime are BPT-SF. 

An important difference between the MEx and CEx di- 
agnostics and our TBT diagnostic is the classification of 
X-ray selected AGNs. While 8% (8/101) of the Juneau 
et al. and 22% (30/126) of the Yan et al. X-ray selected 
L x > 10 42 erg s" 1 AGNs lie in the MEx-SF and CEx- 
SF regimes of their diagnostics, respectively, only 3% 
(3/103) of our X-ray selected AGNs lie within our TBT- 
SF regime (see Section [4]). This may be a result of the 
higher ionization potential of the [Nelll] line (63.45 eV) 
as compared with the [OIII] line (54.94 eV). The higher 
ionization potential appears to foster a more reliable se- 
lection of AGN-dominated galaxies, i.e., the weaker ion- 
izing flux in star-forming galaxies lessens their ability to 
excite [Nelll] as compared to [OIII] . 

We note that our approach does not address the 
issue of classifying galaxies with very low equivalent 
width emission lines, where some of the lines are too 
noisy for traditional line diagno stics to be used. Sec 
ICid Fernandes et all (|2010L 120111 ) for their discussion of 
the WHAN diagram - W^ a versus [Nil] /Ha - and its 
ability to provide a more comprehensive emission line 
classification of galaxies. 

7.1. BPT-comp: AGN-dominated sources 

A critical difference between the iLamareillel ([2010) 
and CEx diagrams and our TBT diagnostic is the lo- 
cation of the BPT-com p (g alaxies that fall bet ween the 
iKauffmann et all 120031 and iKewley et al.l l200l| div isions 
in the BPT diagram). In the lLamareillg (|2010T ) dia- 
gram, ~ 83% of BPT-comp l ie wit hin t he SFG regime. 
In the combined Lamarcillej (|2010f ) and iMarocco et al.l 
(|2011f) diagnostic, ~ 60% of BPT-comp lie within the 
SFG+SFG/comp regimes. In the CEx diagram, ~ 75% 
of BPT-comp lie within the CEx-SF regime. In our 
TBT diagnostic, on the other hand, 65% of the BPT- 
comp lie within the TBT-AGN regime, with only 35% 
in the TBT-SF regim e. Similarly, in the MEx diagnostic 
(| Juneau et all 120111 ). only 17% of BPT-comp lie within 
the MEx-SF regime. 

A number of optical emission-line studies have ar- 
gued that the signal in BPT-comp is dominated by 
star formation activ ity, rather than AGN activity. The 
IKewlev et all (|200ll) upper boundary to the BPT-comp 
regime marks their theoretical prediction for galaxies 
whose contribution from AGN activity to the extreme 
ultraviolet ionizing radiation field just begins to exceed 
50%. According to this work, all sources to the lower- 
left of this boundary have the ir signal domina t ed by star 
formation activity. Similarly, Stasinska et al.l (|2006f ) use 
spectral synthesis modeling to argue that the contribu- 
tion from AGN activity to the emission- line signal in 
BPT- comp is 20% or less. Furthermore, IKewlev et al.l 
(2006) note that BPT -comp lie in the same parameter 
space as HII regions in the [OIII] /[Oil] versus [OI] /Ha 
diagnostic (as well as within the SF-dominated regime 
of the [OIII]/H/3 vs. [SII]/Ha diagnostic). They argue 
that this provides further support for the idea that the 
ionizing radiation field and ionization parameter in BPT- 



comp are dominated by star formation activity. 

Here we argue that our TBT diagnostic's reliance on 
[Nelll], with its higher ionization potential than [OIII], 
[Nil], or [SII], leads to a more reliable identification of 
AGN-dominated sources. In Section[6]we tested whether 
the inclusion of the majority of BPT-comp in our TBT- 
AGN regime was justified. Of the individually X-ray 
detected BPT-comp, 70% are X-ray hard with Lx / £tir 
ratios indicative of dominance by AGN activity (see Fig- 
ure [6|). The stacked signal from X-ray undetected BPT- 
comp is also X-ray hard and exhibits an Lx /-^tir ratio 
> 3 a above the expected range for pure star-forming 
galaxies. This supports our TBT diagnostic inclusion of 
BPT-comp in the TBT-AGN regime and suggests that, 
on average, the X-ray and optical signal in BPT-comp is 
dominated by AGN activity. 

8. SUMMARY 

We have shown that the TBT diagnostic - rest-frame 
(g — z) color versus [NeIII]/[OII] - reliably sepa- 
rates SDSS SF-dominated sources from AGN-dominated 
sources, as classified according to the classic BPT dia- 
gram. Because both [Nelll] and [Oil] are located in the 
blue end of the optical spectrum, we are able to clas- 
sify galaxies using this diagnostic out to z — 1.4. The 
TBT diagnostic provides a significant extension in red- 
shift compared to the BPT diagram (limited in its use 
with optical spectra to z < 0.5) and the more recent 
[OIII] /H/3 -based diagnostic s (limited to z < 0.9 - 'blue 
diagram'. jLamareilld I2010L IMarocco et al.l 120111 : CEx, 
lYan et al.ll20lU MEx. lJuneau et al.ll2011l) . 

We find that the TBT selection of AGNs matches well 
with an X-ray selection of AGNs, with 97% (100/103) of 
our OPTX X-ray selected AGNs lying within the TBT- 
AGN regime. This suggests that the TBT diagnostic 
is more reliable than the BPT diagnostic in identifying 
X-ray selected AGNs, since the BPT diagnostic misiden- 
tifies - 20% of of our OPTX X-ray selected AGNs as 
BPT-SF, i.e., as star-forming galaxies. This may be a 
result of the higher ionization potential of the [Nelll] 
line (63.45 eV) as compared with the [OIII] line (54.94 
eV). The higher ionization potential appears to foster a 
more reliable selection of AGN-dominated galaxies, i.e., 
the weaker ionizing flux in star-forming galaxies lessens 
their ability to excite [Nelll] as compared to [OIII] . 

We perform X-ray and IR stacking analyses of the 
GOODS-N sample of galaxies with accompanying deep 
Chandra imaging to verify the dominance of AGN ac- 
tivity in our TBT-AGN and star formation activity in 
our TBT-SF. We find that the TBT-AGN, on average, 
are X-ray hard with Lx /-^tir > 3 a above the expected 
range for pure star-forming galaxies. Their X-ray hard- 
ness and excess X-ray signal is likely due to obscured or 
weak AGN activity. The TBT-SF, on the other hand, 
are X-ray soft with Lx/L^yr consistent with pure star- 
forming galaxies. 

We perform the same stacking analyses on the BPT 
categories in order to confirm the selection of the ma- 
jority of BPT-comp as TBT-AGN. As expected, the 
BPT-SF are X-ray soft with Lx/Ltir consistent with 
pure star-forming galaxies and the BPT-AGN are X-ray 
hard with Lx/Ltir > 3 a above the expected range 
for pure star-forming galaxies. Of the individually X-ray 
detected BPT-comp, 70% are X-ray hard with Lx / Ltir 
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ratios indicative of dominance by AGN activity. Our 
stacked BPT-comp is significantly detected in both the 
0.5 — 2 keV and 2 — 8 keV bands and is X-ray hard, 
with (r cff ) = 1.0t°:t Furthermore, the stacked BPT- 
comp Lx/Ltyr. is > 3 a above the expected range for 
pure star-forming galaxies. These findings support our 
TBT diagnostic inclusion of BPT-comp in the TBT-AGN 
regime. 

The BPT-comp (individual and stacked) properties 
suggest that, on average, both their X-ray and optical 
signal is dominated by obscured or weak AGN activity. 
This is in contrast to claims by previous optical emission- 
line studies that the signal in BPT-comp is dominated 
by star formation activity. Therefore, we recommend 
that groups carefully consider the impact of excluding 
or including BPT-comp on the interpretation of their re- 
sults. For example, for studies involving determining the 
bolometric contribution from AGN activity or the role 
of AGN activity in galaxy evolution, we advise maximal 
inclusiveness. Since BPT-comp comprise a significant 
percentage of the overall emission-line galaxy population 
(20% of the SDSS DR8 sample), inclusion of the BPT- 
comp would provide a more comprehensive picture of the 
true impact of AGN activity in these studies. 

On the other hand, for metallicity studies of star- 
formi ng galaxies, we adv ise maximal conservativeness 
(e.g.. lTremonti et al.ll2004T ). Emission lines like [OIII] are 
boosted by AGN activity and can masquerade as indica- 
tors of low metallicity, leading to the misinterpretation 
of results. Therefore, in this case, it is optimal to use a 
diagnostic that reliably excludes all AGNs and AGN/SF 
composites. 
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